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We study the three-quark and multi-quark potentials in SU(3) lattice QCD. From 
the accurate calculation for more than 300 different patterns of 3Q systems, the 
static ground-state 3Q potential V3^q'^' is found to be well described by the Coulomb 
plus Y-type linear potential (Y-Ansatz) within 1%-level deviation. As a clear evi- 
dence for Y-Ansatz, Y-type flux-tube formation is actually observed on the lattice 
in maximally-Abelian projected QCD. For about 100 patterns of 3Q systems, we 
perform the accurate calculation for the 1st excited-state 3Q potential V^q ■ by 
diagonalizing the QCD Hamiltonian in the presence of three quarks, and find a 
large gluonic-excitation energy AE^q = VJq' — ^3^*^ of about 1 GeV, which gives 
a physical reason of the success of the quark model. AE3Q is found to be re- 
produced by the "inverse Mercedes Ansatz", which indicates a complicated bulk 
excitation for the gluonic-excitation mode. We study also the tetra-quark and the 
penta-quark potentials in lattice QCD, and find that they are well described by 
the OGE Coulomb plus multi-Y type linear potential, which supports the flux-tube 
picture even for the multi-quarks. Finally, the narrow decay width of penta-quark 
baryons is discussed in terms of the QCD string theory. 



1. Introduction 

Quantum chromodynamics (QCD), the SU(3) gauge theory, was first pro- 
posed by Yoichiro Nambu^ in 1966 as a candidate for the fundamental the- 
ory of the strong interaction, just after the introduction of the "new" quan- 
tum number, "color" In spite of its simple form, QCD creates thousands 
of hadrons and leads to various interesting nonperturbative phenomena 
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such as color confinement^'^'^ and dynamical chiral-symmetry breaking.'' 
Even at present, it is very difficult to deal with QCD analytically due to 
its strong- coupling nature in the infrared region. Instead, lattice QCD has 
been applied as the direct numerical analysis for nonpcrturbativc QCD. 

In 1979, the ffist application^ of lattice QCD Monte Carlo simulations 
was done for the inter-quark potential between a quark and an antiquark us- 
ing the Wilson loop. Since then, the study of the inter-quark force has been 
one of the important issues in lattice QCD.^ Actually, in hadron physics, 
the inter-quark force can be regarded as an elementary quantity to connect 
the "quark world" to the "hadron world" , and plays an important role to 
hadron properties. 

In this paper, we perform the detailed and high-precision analyses for 
the inter-quark forces in the three-quark and the multi-quark systems with 
SU(3) lattice QCD,9'i°'"'i2,i3,i4,i5,i6,i7 ^nd try to extract the proper pic- 
ture of hadrons. 

2. The Three-Quark Potential in Lattice QCD 

In general, the three-body force is regarded as a residual interaction in 
most fields in physics. In QCD, however, the three-body force among three 
quarks is a "primary" force reflecting the SU(3) gauge symmetry. In fact, 
the three-quark (3Q) potential is directly responsible for the structure and 
properties of baryons, similar to the relevant role of the QQ potential for 
meson properties. Furthermore, the 3Q potential is the key quantity to 
clarify the quark confinement in baryons. However, in contrast to the QQ 
potential,^ there were only a few pioneering lattice studies^^ done in 80's 
for the 3Q potential before our study in 1999,^ in spite of its importance in 
hadron physics. 

As for the functional form of the inter-quark potential, we note two 
theoretical arguments at short and long distance limits. 

1. At the short distance, perturbative QCD is applicable, and there- 
fore inter-quark potential is expressed as the sum of the two-body 
OGE Coulomb potential. 

2. At the long distance, the strong-coupling expansion of QCD is plau- 
sible, and it leads to the flux-tube picture. 

Then, we theoretically conjecture the functional form of the inter-quark 
potential as the sum of OGE Coulomb potentials and the linear potential 
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based on the flux-tube picture, 

i<j ' ' 

where Lmin is the minimal value of the total length of the flux-tube link- 
ing the static quarks. Of course, it is highly nontrivial that these simple 
arguments on UV and IR limits of QCD hold for the intermediate region. 
Nevertheless, the QQ potential 1^qq(?') is well described with this form 

^qqW = -^ + W + ^qq- (2) 

For the 3Q system, there appears a junction which connects the three 
flux-tubes from the three quarks, and Y-type flux-tube formation is 
expected. ^'^'^^'^^'^^'^^ Therefore, the (ground-state) 3Q potential is expected 
to be the Coulomb plus Y-type linear potential, i.e., Y-Ansatz, 

^3q" = -^3Q ur^\ + '^SQimin + C3Q, (3) 

where Lmin is the length of the Y-shapcd flux-tubc. 

For more than 300 different patterns of spatially-fixed 3Q systems, we 
calculate the ground-state 3Q potential V^q' from the 3Q Wilson loop W3Q 
using SU(3) lattice QCD^^'^-^'^^'^^ with the standard plaquette action at the 
quenched level on various lattices, i.e., (/3=5.7, 12^ x 24), (/3=5.8, 16^ x 32), 
(/3=6.0, 16^ X 32) and (/? — 6.2, 24^). For the accurate measurement, 
we construct the ground-state-dominant 3Q operator using the smearing 
method. Note that the lattice QCD calculation is completely independent 
of any Ansatz for the potential form. 

To conclude, we find that the static ground-state 3Q potential V^q' 
is well described by the Coulomb plus Y-type linear potential (Y-Ansatz) 
within 1%-levcl deviation. ^"^^^ To demonstrate this, wc show in Fig. 1(a) 
the 3Q confinement potential V^q^^ , i.e., the 3Q potential subtracted by 
the Coulomb part, plotted against the Y-shaped flux- tube length imin- For 
each /3, clear linear correspondence is found between the 3Q confinement 
potential t^q"^ and Lmin, which indicates Y-Ansatz for the 3Q potential. 

Recently, as a clear evidence for Y-Ansatz, Y-type flux-tube formation is 
actually observed in maximally-Abelian (MA) projected lattice QCD from 
the measurement of the action density in the spatially- fixed 3Q system. ^^^^^ 
(See Figs.l (b) and (c).) In this way, together with recent several analytical 
studies, ^^'^^ Y-Ansatz for the static 3Q potential seems to be almost settled. 
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(a) The 3Q confinement potential V^°'^^ 
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i.e., the 3Q potential subtracted by 



Figure 1. 

the Coulomb part, plotted against the total flux-tube length L,nin of Y-Ansatz in the 
lattice unit, (b) & (c) The lattice QCD result for Y-type flux-tube formation in the 
spatially-flxed 3Q system in MA projected QCD. The distance between the junction and 
each quark is about 0.5 fm. 



3. The Gluonic Excitation in the 3Q System 

In this section, we study the excited-state 3Q potential and the gluonic 
excitation in the 3Q system using lattice QCD.^^'^'^ The excited-state 3Q 
potential V^'q ■ is the energy of the excited state in the static 3Q system. 
The energy difference AE^q = V^q — V^q between V^q' and V^q- is 
called as the gluonic-excitation energy and physically means the excitation 
energy of the gluon-field configuration in the static 3Q system. In hadron 
physics, the gluonic excitation is one of the interesting phenomena beyond 
the quark model, and relates to the hybrid hadrons such as qqG and qqqG 
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in the valence picture. 

For about 100 different patterns of 3Q systems, we calculate the excited- 
state potential in SU(3) lattice QCD with 16^ x 32 at /3=5.8 and 6.0 at the 
quenched level by diagonalizing the QCD Haniiltonian in the presence of 
three quarks. In Fig. 2, we show the 1st excited-state 3Q potential V^q' and 
the ground-state potential V^q^'. The gluonic excitation energy AE^q = 
V^Q- — V^q' in the 3Q system is found to be about IGeV in the hadronic 
scale as 0.5fm < Lmin < 1.5fm. Note that the gluonic excitation energy of 
about IGeV is rather large compared with the excitation energies of the 
quark origin. This result predicts that the lowest hybrid baryon qqqG has 
a large mass of about 2 GeV. 




o 

n 0.5 



2 4 6 8 
L„i„ [lattice unit] 



L^,„ [lattice unit] 



Figure 2. The 1st excited-state 3Q potential V3'q'' and the ground-state 3Q potential 
Vg^'^'. The lattice results at /3 = 5.8 and /3 = 6.0 well coincide except for an irrelevant 
overall constant. The gluonic excitation energy AE^q = V^q' — V^q' is about IGeV in 
the hadronic scale as 0.5fm < Lmin < 1.5fm. 



3.1. Inverse Mercedes Ansatz for the Gluonic Excitation in 
3Q Systems 

Next, we investigate the functional form of AE^q = V^q' — V^q' , where 
the Coulomb part is expected to be canceled between V^q' and V3Q ■ . After 
some trials, as shown in Fig. 3, we find that the lattice data of the gluonic 
excitation energy A£'3q = V^q - — V^q' are relatively well reproduced by 
the "inverse Mercedes Ansatz" 

AEsQ = ^ + G, (4) 

3 , 1 
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where Ly denotes the "modified Y-length" defined by the half perimeter 
of the "Mercedes form" as shown in Fig. 3 (a). As for {K, G, we find 
{K ~ 1.43, G ~ 0.77 GeV, ( ~ 0.116 fm) at /3 = 5.8, and {K ~ 1.35, G ~ 
0.85 GeV, C ~ 0.103 fm) at /3 = 6.0. 

The inverse Mercedes Ansatz indicates that the gluonic-excitation mode 
is realized as a complicated bulk excitation of the whole 3Q system. 




2 4 6 8 
Ly [lattice unit] 



5 10 
Ly [lattice unit] 



15 



Figure 3. (a) The Mercedes form for the 3Q system, (b) & (c) Lattice QCD results of 
the gluonic excitation energy Ai?3Q = V^q' — V^q' in the 3Q system plotted against 
the modified Y-length Ltf. The dashed curve denotes the inverse Mercedes Ansatz. 



3.2. Behind the Success of the Quark Model 

Here, we consider the connection between QCD and the quark model in 
terms of the gluonic excitation. ^^'^^'^^'^^ While QCD is described with 
quarks and gluons, the simple quark model successfully describes low-lying 
hadrons even without explicit gluonic modes. In fact, the gluonic excitation 
seems invisible in low-lying hadron spectra, which is rather mysterious. 

On this point, we find the gluonic-excitation energy to be about IGeV 
or more, which is rather large compared with the excitation energies of the 
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quark origin. Therefore, the contribution of gluonic excitations is consid- 
ered to be neghgible and the dominant contribution is brought by quark 
dynamics such as the spin-orbit interaction for low-lying hadrons. Thus, the 
large gluonic-cxcitation energy of about IGeV gives the physical reason for 
the invisible gluonic excitation in low-lying hadrons, which would play the 
key role for the success of the quark model without gluonic modes. ^^'^^'^^'^^ 
In Fig. 4, we present a possible scenario from QCD to the massive quark 
model in terms of color confinement and dynamical chiral-symmetry break- 
ing (DCSB). 



Quantum Chromodynamics 



(Strong) Color Confinement 



Dynamical Chiral Symmetry 
Breaking 



(Tightly Stretched) 
Color Flux- Tube Formation 
with a Large String Tension 
a ~ 1 GeV/fm 



Large Excitation Energy 
of the Flux- Tube Vibration 



Large Constituent Quark Mass 
M„ ~ 300 MeV 



Large Gluonic Excitation 
Energy AE' ~ 1 GeV 



Non-Relativistic Treatment 
on Quark Djaiamics 
in Low-lying Hadrons 



Absence of Gluonic Excitation 
Modes in Low-lying Hadrons 

Massive Quark Model for Low-lying Hadrons 

Figure 4. A possible scenario from QCD to the quark model in terms of color con- 
finement and DCSB. DCSB leads to a large constituent quark mass of about 300 McV, 
which enables the non-rclativistic treatment for quark dynamics approximately. Color 
confinement results in the color flux-tube formation among quarks with a large string 
tension of cr ~ 1 GeV/fm. In the flux-tube picture, the gluonic excitation is described 
as the flux-tube vibration, and its energy is expected to be large in the hadronic scale. 
The large gluonic-excitation energy of about 1 GeV leads to the absence of the gluonic 
mode in low-lying hadrons, which plays the key role to the success of the quark model 
without gluonic excitation modes. 
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4. Tetra-qucirk and Penta-quark Potentials 

In this section, we perform the first study of the multi-quark potentials in 
SU(3) lattice QCD, motivated by recent experimental discoveries of multi- 
quark hadrons, i.e., X(3872) and Z>s(2317) as the candidates of tetra-quark 
(QQ-QQ) mesons, and e+(1540), S — (1862) and 9^(3099) as penta-quark 
(4Q-Q) baryons.^^ As the unusual features of multi-quark hadrons, their 
decay widths are extremely narrow, e.g., r(X(3872)) < 2.3MeV (90 % 
C.L.). For the physical understanding of multi-quark hadrons, theoretical 
analyses are necessary as well as the experimental studies. In particular, 
to clarify the inter-quark force in the multi-quark system based on QCD is 
required for the realistic model calculation of multi-quark hadrons. 

4.1. OGE Coulomb plus Multi-Y Ansatz 

As the theoretical form of the multi-quark potential, we present onc-gluon- 
exchange (OGE) Coulomb plus multi-Y Ansatz^^'^^'^^ based on Eq.(l), i.e., 
the sum of OGE Coulomb potentials and the linear confinement potential 
proportional to the length Lmin of the multi-Y shaped flux-tube. 

On the 4Q potential V4Q, we investigate the QQ-QQ system where two 
quarks locate at (ri, r2) and two antiquarks at (ra, r4) as shown in Fig. 5. 
For the connected 4Q system, the plausible form of V4Q is OGE plus multi- 
Y Ansatz, ^'^ 

K4Q = -^4q{( — + —) + ^( — + — + — + — )} + C^4Qimin + C4Q, (5) 

ri2 r34 2 ri3 ru r23 r24 
while V4Q for the disconnected 4Q system would be approximated by the 
"two-meson" Ansatz as V2QQ = Vq^in^) + VQQ(r24). _ 

On the 5Q potential V5Q, we investigate the QQ-Q-QQ system where 
the two quarks at (ri, r2) and those at (r3, r4) form 3 representation 
of SU(3) color, respectively, and the antiquark locates at r^, as shown in 
Fig. 5. For the 5Q system, OGE Coulomb plus multi-Y Ansatz is expressed 
as V5Q = Vgq'"' + cTsgLmin -I- C5Q with the Coulomb part as 

yCoul ^ _ ^ {(^ + ^) + 1 + J_ + J_ + ^) 

^ ri2 r34' 2Vi5 r25 r35 r45 

1/ 1 1 1 1 XI 

+ -(— + — + — + —)}. 6 

4 ri3 ri4 r23 r24 

We theoretically set (A^q^aiq) and (A^QjCTsq) to be (^3Q,'73q) ~ 
(0.1366, 0.046a"2-) ^i^^ 3Q potential." Note that there is no adjustable 
parameter in the theoretical Ansatze except for an irrelevant constant. 
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Figure 5. (a) A planar tetra-quark (QQ-QQ) configuration, (b) A planar pcnta-quark 
(4Q-Q) configuration, (c) A twisted penta-quark configuration with Q1Q2 ^ Q3Q4. We 
here take di = 0(2 = = <i4 = <i- 



4.2. Multi-quark Wilson loops and multi-quark potentials 

In QCD, the static multi-quark potentials ean be obtained from the corre- 
sponding multi-quark Wilson loops. As shown in Fig. 6, we define the 4Q 
Wilson loop W4Q and the 5Q Wilson loop Wsq^^'^'^'^^ by 



WiQ = -tr(Miii2M2i?i2), 

— ^ ^ahc^a'h'c' ]\~faa' / 



where Li, Ri, M, Mj {i = 1, 2, 3, 4, j = 1, 2) are given by 

Li,Ri,M,Mj = Pexp{ig [ dx^Ai^{x)} G SU(3)^ 

JLi,Ri,M,Mi 



(7) 



(8) 



i.e., Li^Ri^M ^Mj (i ~ 3,4,_7 — 1,2) are line-like variables and Li^Ri (i = 
1,2) are staple-like variables, and Li2,-Ri2 are defined by 



J\o.' o. — _ ^abc ^a' b' c' jj)h' j^cc' R^' ^ = —c^^^c^'^'^'R^^'R'^'^' f9l 
12 2 12' 12 2 1 2 * \ / 

Note that both the 4Q Wilson loop W4Q and the 5Q Wilson loop W5Q are 
gauge invariant. 
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Figure 6. (a) The tetra-quark (4Q) Wilson loop W4Q. (b) The penta-quark (5Q) Wilson 

loop W5Q. 
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We calculate the multi-quark potentials (V4Q, Vsq) from the multi-quark 
Wilson loops {W4Q, W5Q) in SU(3) lattice QCD with /? = 6.0 (i.e., a ~ 
O.lfm) and 16^ x 32 at the quenched level/^'^^'^^ using the smearing method 
to reduce the excited-state components. In this paper, we investigate the 
planar and twisted configurations for the multi-quark system as shown in 
Fig. 5, and show the results for ^1=^2 = ^3 = ^4 = ^ and h\ = h2 = h/2. 

Figure 7 shows the 4Q potential ViQ}"^ For large h, V4Q coincides 
with the energy yc4Q(c', ft-) of the connected 4Q system. For small ft, 
V4Q coincides with the energy V^qq = 2VQQ(ft) of the "two-meson" sys- 
tem composed of two flux-tubes. Thus, we get the relation of V4Q = 
min{V^4Q(rf, /i), 2yQQ(ft,)}, and find the "flip-flop" between the connected 
4Q system and the "two-meson" system around the level-crossing point 
where these two systems are degenerate as T44q(cJ, ft) = 2yQQ(ft). 




^ — • — • — • — • — • — • — • — ^ ^ — • — • — • — • — • — • — • — ^ 

123456789 123456789 
h h 

Figure 7. The 4Q potential V4Q in the lattice unit for the planar 4Q configuration with 
d = 1 (left) and d = 2 (right) as shown in Fig. 5(a). The symbols denote the lattice QCD 
results. The theoretical curves are added for the connected 4Q system (the solid curve) 
and for the "two- meson" system (the dashed curve). 



In Fig. 8, we show the 5Q potential V^q. The lattice data denoted by 

the symbols are found to be well reproduced by the theoretical curve of the 
OGE plus multi-Y Ansatz^^'^^'^"^ with (AaQjCrsq) fixed to be (^3Q,cr3Q) in 
the 3Q potential. 

As a remarkable fact, we find the universality of the string tension and 
the OGE result among QQ, 3Q, 4Q and 5Q systems asiO'"'i2.i3>i4,i5,i6,i7 

<^QQ - <^3Q 0-4Q (75Q, ^^QQ A3Q ~ A4Q =± A5Q. (10) 

This result supports the flux-tube picture on the confinement mechanism 
even for the multi-quark system.^^'^^'^'' 
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Figure 8. The 5Q potential V5Q in the lattice unit for the planar configurations (left) 
and the twisted configurations (right) as shown in Figs. 5(b) and (c). The symbols denote 
the lattice QCD results. We add the theoretical curve of the OGE plus multi-Y Ansatz 
with (Asq.ctsq) fixed to be (A3q,(T3q). 

5. QCD String Theory for the Penta-Quark Decay 

Our lattice QCD studies on the various inter-quark potentials indicate the 
flux-tube picture for hadrons, which is idealized as the QCD string model. 
In this section, we consider penta-quark dynamics, especially for its ex- 
tremely narrow width, in terms of the QCD string theory. 

The ordinary string theory mainly describes open and closed strings cor- 
responding to mesons and glueballs, and has only two types of the reaction 
process as shown in Fig. 9: 

1. The string breaking (or fusion) process. 

2. The string recombination process. 




Figure 9. The reaction process in the ordinary string theory: the string breaking (or 
fusion) process (left) and the string recombination process (right). 



February 1, 2008 4:48 WSPC/Trim Size: 9m x 6m for Proceedings 



Conf6 



12 

On the other hand, the QCD strmg theory describes also baryons and 
antibaryons as the Y-shaped flux-tube, which is different from the ordinary 
string theory. Note that the appearance of the Y-type junction is peculiar 
to the QCD string theory. Accordingly, the QCD string theory includes the 
third reaction process as shown in Fig. 10: 

3. The junction (J) and anti-junction (J) par creation (or annihilation) 
process. 




Figure 10. The junction (J) and anti-junetion (J) par creation (or annihilation) process 
pecuUar to the QCD string theory. 

Through this J- J pair creation process, the baryon and anti-baryon pair 
creation can be described. 

As a remarkable fact in the QCD string theory, the decay process (or 
the creation process) of penta-quark baryons inevitably accompanies the J- 
J creation^^ as shown in Fig.ll. Here, the intermediate state is considered 
as a gluonic-excited state, since it clearly corresponds to a non-quark-origin 
excitation. 

Our lattice QCD study indicates that such a gluonic-excited state is 
a highly-excited state with the excitation energy above IGeV. Then, in 
the QCD string theory, the decay process of the penta-quark baryon near 
the threshold can be regarded as a quantum tunneling, and therefore the 
penta-quark decay is expected to be strongly suppressed. This leads to a 
very small decay width of penta-quark baryons. 




Figure 11. A decay process of the penta-quark baryon in the QCD string theory. The 
penta-quark decay process inevitably accompanies the J-J creation, which is a kind of 
the gluonic excitation. 



Now, we try to estimate the decay width of penta-quark baryons in 
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the QCD string theory. In the quantum tunnchng as shown in Fig.ll, 
the barrier height corresponds to the gluonic excitation energy A_E of the 
intermediate state, and can be estimated as ~ IGeV. The time scale 
T for the tunnehng precess is expected to be the hadronic scale as T = 

0. 5 ~ Ifm, since T cannot be smaller than the spatial size of the reaction 
area due to the causality. Then, the suppression factor for the penta-quark 
decay can be roughly estimated as 

I exp(-A£T)|2 ~ |e-lGeVx(0.5^1)fm|2 ^ ^ ^Q-i (11) 

Note that this suppression factor | exp(— A£^T)|^ appears in the decay pro- 
cess of the penta-quark baryons for both positive-parity and negative-parity 
states. 

For the decay of e+(1540) into N and K, the Q-value QisQ = M(e+)- 
M(N) - M(K) ~ (1540 - 940 - 500)MeV ~ lOOMeV. In ordinary sense, 
the decay width is expected to be controlled by Fhadron — Q — lOOMeV. 
Considering the extra suppression factor of | exp(— Ai^T)^, we get a rough 
order estimate for the decay width of e+(1540) as r[6+(1540)] ~ Thadron x 
|exp(-A£;r)|2 ~ 1 ~ IQ-^MeV. 
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